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Possibility of gravitational repulsion in General Relativity is discussed and astro- 
nomical data in favor of cosmological acceleration are described. The problem of 
vacuum energy is emphasized and possible ways of its solution are indicated. The 
main attention is payed to adjustment mechanism which in principle could compen- 
sate originally huge vacuum energy down to cosmologically acceptable value and to 
solve the coincidence problem of a close magnitudes of the non-compensated rem- 
nants of vacuum energy and the energy density of the universe at the present time. 
Finally possible modifications of gravity at large scales which could induce acceler- 
ated cosmological expansion are considered. 



I. INTRODUCTION 

Astronomical data accumulated during the last two or three decades unambiguously 
proved that the universe expands with acceleration, i.e. the cosmological expansion velocity 
rises with time. To appreciate the surprising features of this result let us consider a simple, 
though not precise, analogy. A stone thrown vertically up from the Earth surface would 
either stop at some moment and return back to the Earth or, with sufficiently large initial 
speed, would fly away forever to infinity. There is a special intermediate case when the 
velocity of the stone at infinity tends to zero. In all the cases the acceleration is directed 
towards the Earth, i.e. the velocity during "expansion" drops down. It was believed until 
recently that the cosmological expansion proceed qualitatively in the same way. The speed 
of the expansion decreases with time and either the expansion will stop and the universe 
will collapse back to a hot singular state, or the expansion will last forever. As in the 
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case of the stone considered above the velocity of the expansion always drops down with 
time. The intermediate situation of zero velocity at infinity (i.e. expansion stops at infinity) 
corresponds to spatially flat universe with 3D Euclidean geometry. In this old picture 
the spatial geometry and the ultimate destiny of the universe are rigidly connected. If the 
universe is spatially closed (3D sphere), the expansion turns into contraction. If the universe 
is open (3D hyperboloid) the expansion never stops. 

It is understood now that this picture is not true. First, as is mentioned above, the 
universe today expands with acceleration but most of its previous history (in terms of the 
cosmological red-shift but not time, see below) it expanded with the normal deceleration. It 
is similar to a stone thrown up from the Earth surface which first moved with a decreasing 
velocity but after a while it started to accelerate (as if it has a rocket accelerator), gained 
speed, and never came back. It is practically established that at the very beginning the 
universe also expanded with acceleration. It is the famous inflationary stage when the initial 
cosmological push was operated for a short while. It is analogous to the initial acceleration 
of the stone in the considered example. After that the subsequent motion both of the stone 
and of the universe was simply the inertial one. 

As we understand, most probably the initial inflationary push and the cosmological ex- 
pansion today are both created by an antigravitating state of matter, though some other 
mechanisms are not excluded. The antigravity at the beginning is necessary for our existence 
because it was a source of expansion which created large and suitable for life universe but 
it seems unnecessary now, or we do not understand its necessity. 

The paper is organized as follows. In the next section the basic cosmological equations are 
derived in a simplified way and different expansion regimes are considered. In sec. 3 main 
cosmological parameters are introduced and the data in favor of the accelerated expansion 
are presented. In sec. 4 the problem of vacuum energy is described and some suggestions for 
its solution are discussed. In sec. 5 modified gravity models leading to accelerated expansion 
are considered. 

II. BASIC EQUATIONS 

The distribution of matter in the universe is assumed to be homogeneous and isotropic, 
at least at an early stage, as indicated by the isotropy of the Cosmic Microwave Background 
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radiation (CMB), and even now at very large scales, larger than 100 Mpc. Correspondingly 
the metric can be chosen in the Friedman- Robert son- Walker (FRW) form: 

ds 2 = dt 2 - a 2 {t) \f{r)dr 2 + r 2 dVt] , (1) 

where f(r) describes three-dimensional (3D) space of constant curvature, f(r) = 1/(1 — kr 2 ) 
with k — +1 corresponding to 3D hyperboloid, k = — 1 corresponding to 3D sphere, and 
k = corresponding to flat 3D space. 

The evolution of scale factor a(t), i.e. the expansion law, is determined by the Friedmann 
equations, which follow from the general General Relativity (GR) ones: 

Rill/ ~^.9uvR = 5~ Tn V (2) 

2 m Pl 

for FRW anzats (TO . The term in the r.h.s. of this equations is the energy- momentum tensor 
of matter, which is the source of gravity in GR. This equation has a very important property 
that its left hand side (l.h.s.) is covariantly conserved: 

D, (K ~ \^R) = , (3) 

where is the covariant derivative in the corresponding Riemann space-time. This identity 
is the so called Bianchi identity valid in any metric theory. The conservation of the l.h.s. 
implies the conservation of the r.h.s.: 

DPS = 0. (4) 

On the other hand this law of covariant conservation of the energy-momentum tensor is 
fulfilled due to invariance of the theory with respect to general coordinate transformation, 
i.e. general covariance. In this sense the theory is self-consistent. 

The derivation of the Friedmann equations from Einstein equations (j2J) is straightforward 
but quite tedious. So we proceed in a simplified, non-rigorous way which not only allows 
to avoid complicated algebra but gives an intuitive understanding of the equations. Let us 
consider a small sphere in the cosmological matter with radius a and a test body on the 
surface of this sphere. The matter outside this sphere does not have any gravitational action 
on the test body, so the law of the energy conservation for such test body can be written as 

E kin + U = v 2 /2 -G N M/a = const , (5) 
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where Gn = l/ m p« * s the gravitational coupling constant, M = 4na 3 p/3 is the total mass 
inside radius a, and p is the energy density of matter inside the sphere. 

Introducing the Hubble parameter H = a/ a we rewrite this equation in the canonical 
form: 



H 2 _ ia\ 2 _ 8npG N k ^ 
a) 3 a? 



From this definition and the form of metric (CQ) follows the Hubble expansion law: 

v = Hd, (7) 

where v is the velocity of two non-interacting objects at distance d. Note that objects in 
bound systems e.g. (stars in galaxies)) do not flow away from each other. The Hubble 
expansion law is simply motion by inertia of previously accelerated and now non-interacting 
(except for gravity) objects. 

One more equation follows from the energy balance in volume element V of the medium: 
dE = —PdV, where P is the pressure density and E = pV, so that dE = Vdp + 3(da/a)Vp. 
Hence: 

p + 3H(p + P) = 0. (8) 

This equation is the condition of the covariant conservation of the energy-momentum tensor 
01]) for a special case of FRW metric. 

From eqs. (E]) and (jSJ) follows the expression for the acceleration of the test body: 

d ATT Gat , „. 

- = —^(P + 3F). (9) 

Let us stress that not only energy but also pressure gravitate and induce gravitational 
acceleration. Here lays a fundamental difference between the General Relativity and the 
non-relativistic Newtonian theory of gravity. 

It looks quite surprising that our non-relativistic analysis of the behavior of the test body 
in seemingly Newtonian gravitational theory have lead us to GR equation (Q. The essential 
input, where GR is involved, is the assumption that the source of gravity is not mass but 
energy. Together with the law of the energy balance/ conservation (jBJ) it leads to participation 
of pressure in creation of the gravitational force and this is of fundamental importance for 
existence of life. The energy density of matter must be positive in any known reasonable 
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non-pathological theory. So if only energy density gravitated, gravitational repulsion would 
be impossible and the universe would not expand to a possible-for-life state. However, 
pressure may be negative and if P < — p/3, acceleration a would be positive and it would 
result in cosmological antigravity and induce cosmological expansion. Life is impossible with 
Newtonian gravity. 

Existence of gravitational repulsion in cosmology does not mean that it may be possible 
to construct an anti-gravitating space ship. Antigravity is allowed only for infinite size 
objects, such as the universe as a whole, domain walls, or cosmic strings. One can see that 
finite objects can create only gravitational attraction analyzing the Schwarzschild solution, 
according to which the gravitational field created by a compact object is proportional to a 
positive definite integral of the energy density multiplied by a certain positive function which 
includes the gravitational mass defect. The mass, which creates the gravitational field, 



is positive for any pressure. 

Another way to check that the pressure inside a compact source does not gravitate is to 
integrate by parts the expression: 



To complete the set of the equations describing the cosmological evolution we need to 
know the equation of motion of a test particle in the FRW space-time. If 3D space is flat, 
the variation of the particle momentum is simply determined by the Doppler effect. Indeed, 
let us consider a particle moving from point x\ to 22 separated by distance dx — X2 — X\. 
According to the Hubble law fl7D, the relative velocity of these two points is U — Hdx. 
Correspondingly the Doppler shift of the particle momentum is dp = —UE = —HEdx and 
so: 



In other words, the momentum of a free particle drops down as the inverse scale factor. At 
this stage it is convenient to introduce the cosmological redshift: 




(10) 




(11) 



p = —HEx = —Hp . 



(12) 



z = a(tu)/a(t) — 1 



(13) 
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So the solution of eq. (fT2l) looks as: p ~ 1/a ~ l/(z + 1). 

Eq. (fT2l) is the geodesic equation in 3D flat FRW metric. The general form of the geodesic 
equation is: 

dV a 

_ = -r^vv , (14) 

where = dx a /ds and r° are the Christoffel symbols in FRW metric. If the 3D space 
is not flat, there appeared some additional terms proportional to spatial curvature, which 
usually are negligibly small. 

Out of three equations ([8]), and ([9]) only two are independent but they contain three 
unknown functions of time: p, P, and a. Hence for the description of the cosmological 
evolution one more equation is necessary. Normally to this end the equation of state is 
used defining pressure density in terms of energy density: P = -P(p). In many practically 
interesting cases the linear equation of state is sufficient: 

P = wp (15) 

where w is (normally) a constant parameter, different for different forms of matter. 

Now we discuss three equations of state which were realized in cosmology at different 
periods of the universe evolution. 

For nonrelativistic matter P <C p . Hence in a good approximation one can take w — 0. 
So from eq. (IE]) follows pNR = SHp^ji and thus: 

p~l/a 3 . (16) 

The result is physically evident: the number density of the particles drops down as volume 
(as it is sometimes said that the number density is conserved in comoving volume, i.e. in 
the volume expanded together with the universe), while the energy of each nonrelativistic 
particle is equal to its mass and remains constant. 

The expansion law can be found from equation fl6]) and in the particular case of k = 
has a very simple form: 

a NR ~ t 2/3 - (17) 

For relativistic matter P = p/3. This expression is intuitively evident because for rela- 
tivists particle momentum is equal to its energy and the pressure is obtained by averaging 
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of the momentum flux over three spatial directions. Correspondingly, according to eq. (JHJ) 
the energy density evolves as 

Prel^l/a 4 . (18) 

An extra factor 1/a in comparison with nonrelativistic matter appears because of the en- 
ergy/momentum redshift of individual relativistic particles. Correspondingly in relativistic 
expansion regime a ~ t 1 ^ 2 . 

Another very interesting and rather unusual equation of state is the so called vacuum(- 
like) one, when w = —1, which is (approximately) realized in the universe at the present 
time. Its source can be in particular the vacuum energy-momentum tensor. Since vacuum is 
supposed to be the same independently of the reference frame, its energy-momentum tensor 
must have the same form in any reference frame too and thus it must be proportional to the 
only "invariant" second rank tensor, i.e. to metric tensor, g^: 

P^iv Pvac 9fiv and Pvac Pvac • (^^) 

Since p = —3H(p + P) = 0, vacuum energy density remains constant in the course of the 
cosmological expansion: 

Pvac = const . (20) 

With constant energy density the Hubble parameter is constant too and the expansion is 
exponential, a ~ exp(H vac t). 

Such strange behavior of the energy density which does not decrease, while the volume 
rises exponentially, leads to a very important and striking phenomenon. Namely if in some 
piece of the "primordial" universe, the equation of state would be the vacuum- like one, 
( I19p . this presumably microscopically small volume would quickly expand becoming larger 
than the present-day universe and all the matter in our universe would be created from the 
stretched out microscopically small amount of matter in the original volume. So all our 
visible universe could originate from microscopically small volume with negligible amount 
of matter by exponential expansion with constant p. 

Note that if k > and p ~ l/a n with n=3,4, the cosmological expansion will change to 
contraction. If k < 0, the expansion will be eternal for any equation of state. However, if 
p > k/a 2 , as e.g. p vac = const, the expansion will last forever for any k. 
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As we have seen, the set of equations (jSJ), ©, and (jHJ) does not have stationary solutions. 
It was to great disappointment of Einstein who wanted to have static universe despite the 
problem of thermal death of the universe and the Olberts paradox. So Friedman [l] , who first 
found the non-stationary realistic cosmological solutions predicted the universe expansion, 



later discovered by Le Maitre [2] and Hubble |3|]. Thus the Hubble law would be more 
appropriate to name the Friedman-LeMaitre- Hubble law (FMH law). After the works of 



Robertson and Walker [4] , who studied general properties of the metric associated with the 
Friedman solution, it got the name FRW metric. 

Note in conclusion to this section that if w < —1, the energy density rises with time and 
such behavior leads to the so called phantom singularity: 

, i v2(M-l)/3 , 1 N 2(|iu|-1) 

a «~G^) ■"~(^) ' (21) 

Thus in finite time both the scale factor and the energy density would become infinitely 
large and not only galaxies, stars, and planets, but atoms and even elementary particles 
would be turn apart. The existing astronomical data indicate to the possibility that w < — 1 
but it does not necessarily means that the universe will end in singular state because w 
may change with time and return to w > —1 before the singularity is reached. For example 
w(t) < —1 may be realized by a scalar field with some unusual form of the kinetic term, but 
in the course of the field evolution (w + 1) may become non-negative. 

The energy-momentum tensor of a scalar field with the normal kinetic term (f) has the 
form: 

T^d = d^d v <P - [(«90) 2 - U {<!>)] (22) 

Correspondingly for a homogeneous scalar field </> = 4>(t) the energy and pressure density 
are: 

, = £+w P = t -uw (23) 

and the condition \P\ < p is always fulfilled, if pressure is negative. Clearly for a slowly 
varying the vacuum-like equation of state P ~ — p is realized. In the limit of very fast 
field variation when the potential can be neglected the most rigid equation of state P = p 
becomes valid. 

Note that the equation of state in the form P = P{p) does not always exist. The 
relation between P and p may be non-local in time, though one can always define parameter 
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w(t) = P(t)/p(t). Despite an absence of the equation of state, the set of cosmological 

equation can be still complete because the necessary missing equation is the equation of 
motion for <p\ 

D 2 <f) + U'(<f)) = 0. (24) 

III. UNIVERSE TODAY 

The rate of the universe expansion, see eq. (JTj), is determined by the present day value 
of the Hubble parameter: 

H = 100 h km/sec/Mps , (25) 

where h = 0.73 ± 0.05 and H' 1 = 9.8 Gji/h « 13.4 Gyr. 

The cosmological energy densities of different forms of matter, p a are expressed through 
dimensionless parameters: 

= Pa/Pc, (26) 

where p c is the so called critical energy density, i.e. energy density in spatially flat universe 
corresponding to k = 0. From equation (jHJ) it follows: 

p c = ^H3l = 1.88 • 10~ 29 /i 2 ^ = 10.5 h 2 ^, « 10~ 47 h 2 GeV 4 . (27) 
87T cm d cm-* 

It corresponds to 10 would-be protons per m 3 , but in reality the density of protons is much 
smaller; the dominant form of matter in the universe is not the usual baryonic one but 
something invisible called dark matter and dark energy. 

Matter inventory (for a review see e.g. ref. 5[): 
The total cosmological energy density is quite close to the critical one: 

n to t = $>a/p c = l±0.02. (28) 

a 

This result is obtained from the analysis of the angular fluctuations of CMBR, in particular, 
from the the position of the first acoustic peak peak and from the study of the large scale 
structure (LSS) of the universe. 

The usual baryonic matter contributes a minor fraction to the cosmological energy density, 
Qb — 0.044 ± 0.004. It is measured by comparison of the relative heights of CMBR peaks 
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and by the analysis of big bang nucleosynthesis (BBN). Qualitatively one can conclude that 
the fraction of baryons is small because in the baryon dominated universe the structure 
formation could start only rather late, after recombination at z ~ 10 3 , but in this case there 
is too little time for density perturbations to evolve up to the observed values, Sp/p > 1. 

Thus 95% of matter in the universe is something unknown, dark. This dark staff can be 
separated into two quite different components: dark matter (DM) and dark energy (DE). 
Dark matter is supposed to consist of some unknown objects (new stable elementary par- 
ticles, primordial black holes, topological or non-topological solitons, etc). The fraction of 
dark matter is five times larger than fraction of baryons: Qdm ~ 0.22 ± 0.04. This number 
was obtained from galactic rotation curves, gravitational lensing, equilibrium of hot gas in 
rich galactic clusters, cluster evolution, and baryon acoustic oscillations. The equation of 
state of dark matter is simple non-relativistic one with w — 0. 

The remaining three fourth of the cosmological energy density consists of mysterious 
dark energy with Qde ~ 0.75 and negative parameter of the equation of state, w « — 1. 
According to eq. ([9]), such state of matter induces accelerated expansion. The cosmological 
fraction of DE was found from the dimming of high-z supernovae, CMBR, LSS, and from 
the universe age. 

We can calculate the universe age in terms of the present day values of the cosmological 
parameters integrating the equation: 

r n 1/2 

a= [8tc pG N a 2 /3- k\ . (29) 

As a result one obtains 

tu = 77 / / > (3°) 



where Q m and Q r are the fractions of non-relativistic and relativistic matter respectively, 
and fl v is the fraction of vacuum energy with P = —p. 

According to the ages of old stellar globular clusters and nuclear chronology the universe 
age lays in the interval: 

t v = 12- 15Gyr, (31) 

which well agrees with the presented above data on H and fl a but disagrees with decelerated 
universe without dark energy. 
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Let us stress that different methods of determination of the cosmological parameters are 
absolutely independent. They are obtained by different types of astronomical observations 
and rely on different physical phenomena in the universe. The latter excludes possible 
interpretation errors. So phenomenologically it is established without any doubts that the 
universe expands with acceleration and the source of this accelerated expansion makes 75% 
of the total cosmological energy density but it is unknown what exactly it is. There are two 
main competing possibilities: a new light or massless (scalar) field or modification of gravity 
at cosmologically large scales. The first one includes in particular some not completely 
compensated remnant of vacuum energy carried by the compensating field. In the case of 
modified gravity one must take care not only of the dimming of the high ref-shift supernovae 
but also of the large scale structure formation, CMB fluctuations, universe age, etc. 



IV. PROBLEM OF VACUUM ENERGY 



The problem of vacuum energy is probably the most striking problem of contemporary 
fundamental physics. This is a unique example when theoretical expectations differ from the 
observations by 50 -100 orders of magnitude. The beginning of the story is traces back to 
almost century ago [6] when Einstein introduced into GR equations ([2]) an additional term 
proportional to metric tensor g^: 

R^u - ^g^R - kg^u = 8ttG n . (32) 

Coefficient A must be constant to satisfy the constraints of general covariance and energy- 
momentum conservation, see eqs. fl3]) and (jlj). By this reason Lambda-term is often called 
cosmological constant. 

Evidently A-term is equivalent to vacuum energy f Tl9|) . though there are still some state- 
ments in the literature that A is a geometrical quantity, while T^ ac ) has completely different 
nature and so they should be treated differently. However, there is no way to distinguish 
them observationally and thus they are absolutely the same from physical point of view. 

As follows from quantum field theory, vacuum energy of any quantum field is infinitely 
large. Indeed quantum field can be understood as a collection of quantum oscillators with 
all possible frequencies each having the ground state energy equal to u/2. Integration over 
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all frequencies for a bosonic field gives: 

Pvlc = (Kb) vac = [ 7^ ^(4°fe + b kbl) V ac = [ 7^ U k = OO* . (33) 



(2vr) 3 2 ™ " ' K klvac J (2tt) 



.(/) 

rvac 



{U f ) vac = J y (4«fc - fefcftfc)^ = / 7^4i w fc = -oo 4 . (34) 



Here = \/k 2 + m 2 and m is the mass of the field. 

On the other hand, fermionic vacuum fluctuations also have infinitely large vacuum energy 
but of the opposite sign: 

d 3 k u k t f /• rf 3 fc 

(2tt) 3 2 l%afc Wmc ~ 7 (2vr) : 

It was noticed by Zeldovich [7| that in the world with an equal number of bosonic and 
fermionic species having equal masses, at least pairwise, energy of vacuum fluctuations 
vanishes. Later it was suggested that there may indeed exist symmetry between bosons and 
fermions called supersymmetry (SUSY) |8[, which demands an equal number of bosonic and 
fermionic degrees of freedom. If so, then the quartically divergent contribution to vacuum 
energy vanishes. Clearly supersymmetry is not exact because there are no supersymmetric 
partners of the observed particles with the same masses. Presumably they are much heavier, 
at least by msusY ~ 10 3 — 10 4 GeV. In this case the sum of bosonic ( |33l and fermionic (|34l) 
contributions would be "only" quadratically infinite, p^} c + p\/J c ~ m %usY x °° 2 - With a 
particular mechanism of SUSY breaking the inifinities may be cancelled but the net result 
must be non-zero and of the order of Tng USY , which is larger than the observed value of 
pvac by at least 56 orders of magnitude. In the version of supersymmetry which includes 
gravity (supergravity) the net result for vacuum energy may vanish in the broken case, 
but at the expense of a fantastic fine-tuning. The natural value of p vac in supergravity is 
about Wpj ~ 10 76 GeV 4 , so the fine-tuning of the parameters must be achieved with the 
unbelievable precision 10~ 123 . 

Let us note also that vacuum (or vacuum-like) energy undergoes several huge jumps in 
the course of the cosmological evolution. At the very beginning during inflationary epoch 
the vacuum-like energy of the inflaton field could be as large as 10 56 GeV 4 . Though it was 
not real vacuum energy but still the inflaton was anti-gravitating similar to the present 
today dark energy. Later in the process of the cosmological cooling down a sequence of 
phase transitions occurred in the primeval plasma. If at the end of inflation the plasma was 
heated up to the temperatures above the grand unification scale, T ~ uigut ~ 10 15 GeV, 
the GUT symmetry would be unbroken and broke downs later at smaller T. The phase 
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transition from symmetric to broken symmetry state is accompanied by the change of the 
ground state energy- momentum tensor approximately equal to AT^ U w ^Gi/rfiV- Similar 
phase transitions took place at the electroweak and QCD scales with change of vacuum 
energy 10 8 GeV 4 and 10~ 2 GeV 4 respectively. There might also be phase transitions related 
to supersymmetry breaking. 

These simple estimates show that vacuum energy is naturally expected to be either infinite 
or almost infinite. The attitude of the community, at least of that part of the community 
which appreciated this problem was quite strange, namely it was implicitly taken that oo = 0. 
Such approach reminds well known quotation from Feynman about radiative corrections in 
quantum electrodynamics: "Corrections are infinite but small" . 

The attitude of the astronomers and cosmologists to Lambda term was (and still is) 
strongly polarized. The dominant part of the establishment believed that Lambda is identi- 
cally zero, while some other thought that introduction of Lambda term was a very important 
generalization of GR (9)]. Among the last was Lemaitre who said that it was the greatest 
discovery, worth alone to make Einstein's name famous. A strong antagonist of Lambda was 



Gamow, who wrote in his autobiography book "My world line" |10|] that A raises its nasty 
head again" after indication at the beginning of the 60s that quasar are accumulated near 
z = 2, which later found "Lambda-less" explanation. 

Now the gravity of the vacuum energy problem placed it into the central position in fun- 
damental physics. Probably the most serious argument in favor that something mysterious 
happens in vacuum comes from quantum chromodynamics (QCD). According to this well 
established theory which beautifully passed all experimental tests, u and d quarks are very 
light. Their masses are about 5 MeV. Proton is known to be a bound state of these three 
quarks, p = (uud). So the proton mass should be 15 MeV minus binding energy, instead of 
938 MeV. The solution of the problem suggested by QCD is that vacuum is not empty but 
filled with quark jll] and gluon 12] condensates: 



m ± o, 

(G^Gn + 0, (35) 

having negative vacuum energy: 

P^ D) « ("0.01 GeV 4 ) » (-10 45 Pc ) • (36) 
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Vacuum condensate is destroyed around quarks and the proton mass becomes 

m p = 2m u + m d - p^ D) l 3 p ~ 1 GeV , (37) 

where l p ~ (O.lGeV -1 ) is the proton size. 

The value of the vacuum energy of the quark and gluon condensates (1361) is practically 
established by experiment. To adjust the total vacuum energy down to the observed magni- 
tude, ~ 10~ 47 GeV 4 , there must exist another contribution to vacuum energy of the opposite 
sign and equal to the QCD one with precision of one part to 10 45 . This new field cannot 
have any noticeable interactions with quarks and gluons, otherwise it would be observed 
in direct experiment, but still it must have very same vacuum energy. This is one of the 
greatest mysteries of Nature. 

At this place it is proper to make an intermediate summary. 

1. There are known and fantastically huge contributions to p vac but a mechanism of their 
compensation down to (almost) zero remains mysterious. 

2. The observed today vacuum energy is very close the energy density of the cosmological 
matter, p vac ~ 3(pb + Pdm), though they evolve in a different way during cosmological 
history, i.e. p ma tter ~ l/^ 2 and p vac = const. Why their near equality happened to occure at 
the present epoch? 

3. What is the nature of anti-gravitating dark energy? It seems to have the equation of 
state consistent with w = —1. Is it simply vacuum energy or something more interesting? 

Mostly only problems 2 and 3 are addressed theoretically (phenomenologically) either by 
(infrared) modification of gravity or by a new (scalar) field (quintessence) leading to the 
accelerated expansion. However evidently all three problems are strongly coupled and can 
be solved only after adjustment of "infinitely" large p vac down to tiny p c is understood. 



A. Some suggestions for the solution of the vacuum energy problem 

There are several suggestions in the literature for possible resolution of the vacuum energy 



problem, for a review see ref . [13| . Below we briefly discuss some of them and dwell in some 
detail on the dynamical adjustment mechanism. The short and incomplete list of different 
option looks as follows: 

1. Subtraction constant. If we encounter exact vacuum energy problem, i.e. w — — 1, 
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a trivial resolution of the difficulty could be a very precise choice of the zero level of the 
energy density or, in other words, the choice of the subtraction constant whose value must be 
taken such that it would very accurately compensate all the contributions to vacuum energy 
from vacuum fluctuations, phase transitions, etc, being taken with opposite sign to them, 
but with small "imperfection" exactly equal to the observed value of p vac . It is impossible 
to exclude such "solution" but it is extremely unnatural. 

If w ytz — 1 ; dark energy is not constant but evolves during cosmological history and simple 
subtraction cannot eliminate it. 



2. Anthropic principle 14|. In some sense it is similar to the solution with help of the 



subtraction constant. Unnaturalness of the latter is lifted by the assumption that there is a 
huge number of different universes with all possible sets of the subtraction constants. The 
conditions suitable for life exist only in universes with rather narrow span of values of p vac . 
Large and positive p vac would produce too fast cosmological expansion, so the structure for- 
mation would be strongly inhibited. The opposite case of negative and large by the absolute 
value p vac would force the universe to collapse back to hot singularity before it becomes suit- 
able for life. An existence of many universes with different physical laws and possibly with 
different subtraction constants in connection with the vcuum energy problem was suggested 
in ref. 15]. Multivacuum states can be realized in brane landscape theories [lij with about 



10 1000 different vacuum states. Large number of vacuum states naturally appear in infla- 
tionary scenarios jl?j |. especially in chaotic inflation 18]. Realistic way of compactification 
in multidimensional brane theories leading to huge number of vacuum states was suggested 
in ref. 19 ] . 

For many people (including the author) the anthropic solution looks unsatisfactory 
because the theory is impossible to falsify. It reminds the situation with the problems 
of the Friedmann cosmology prior to the suggestion of the inflationary paradigm. Before 
inflation the anthropic principle was believed to be the only feasible way to understand how 
or why our universe was created. Inflation not only resolved all the problems in a simple 
and natural way but also predicted some observational effects, in particular, the shape of 



the perturbation spectrum 



which is confirmed by the data, and cosmic background of 



relic gravitational waves 21], which may be a final proof of inflation. 



3. Dynamical adjustment. The idea is similar to axionic solution of the strong CP 
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problem. A new light field is introduced whose back-reaction cancels out the impact of the 
source which created this field. This principle is known in chemistry and physics from the 
XIX century and has the name "the Le Chatelier principle". 

The original [22| and many subsequent suggestions [13[ were based on a scalar field but 
higher spin fields: vector 24] or tensor (25], are also possible. The main idea is that a new 
field $ is coupled to gravity in such a way that vacuum energy through its gravitational 
action leads to formation of the condensate of The energy density, p$, compensates the 
original vacuum energy. If p vac > 0, then p$ should not be positive definite but the related 
instability should not be catastrophic i.e. the energy of $ does not quickly drop down to 
negative infinity. 

Independently on a concrete realization of the idea the generic predictions are similar 
and quite attractive. First, when the contribution of $ on the evolution of the scale factor 
is taken into account, the exponential cosmological expansion turns into the power law one. 
Second, the compensation of the vacuum energy is not complete but it is compensated 
only down to the terms of the order of p c {t). Third, the non-compensated energy may 
have quite an unusual equation of state. Such mechanism not only solves the problem of 
compensation of the original vacuum energy, which is estimated to have a typical value of the 
particle physics scale, down to a cosmologically small value and also explains the so called 
coincidence problem, i.e. the close proximity of the observed p vac to the time dependent pJf . 
This is exactly what is observed. In this sense the dark energy was predicted in 1982 22]. 
Unfortunately, despite the numerous attempts no realistic model was found starting from 
1982 till now. 

Probably the first person who suggested that vacuum energy may be of the order of 
the time dependent cosmological energy density was Matvey Bronshtein 23(. However, 
taken literally, the models with A = A(t) are not innocent because the transversality of the 
Einstein equations demands A = const, see eq. ( !20l . To achieve such transversality, or what 
is the same, the conservation of (jlj), some new light or massless fields are necessary. 
Another possibility demands serious modifications of the theory, e.g. non-metric theories. 
This was a reason why the Bronstein suggestion was strongly criticized by Landau. However, 
an approximate relation P —p can be achieved with light scalar field which might have 
p = p(t), see eq. fl23l) . 

The first model of dynamical adjustment suggested in ref. 22| was based on non- minimally 
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coupled scalar field satisfying the equation of motion: 

4> + 3H<j) + U'((f),R) = 0. (38) 

The non-minimal coupling to gravity was taken in the simplest form U = £,R<p 2 /2. It is 
easy to see that the solutions of this equation are unstable if C.R < because in De Sitter 
background this term behaves as tachyon with negative effective mass squared, < 0. 

One can find that in the initially de Sitter state with the exponentially rising scale factor, 
a(t) ~ exp(H v t), cf) also rises exponentially but when the back reaction of the rising <ft on 
the cosmological expansion is taken into account, the exponential rise turns into the power 
law one. Asymptotically: 

a(t)~t^, (39) 

where /3 is a constant. In other words the initial De Sitter space-time becomes the Friedmann 
one, despite non-zero vacuum energy. Still the vacuum energy is not compensated because 
the energy-momentum tensor of is not proportional to the metric tensor 

T^(0)^A^ (40) 

but the change of the regime is achieved due to the weakening of the gravitational coupling: 

G N ~ 1/t 2 . (41) 

If such rise of Mp\ took place in the early universe and was stabilized at some later 
but still an early stage, this mechanism might explain hierarchy between electro-weak and 



gravitational mass scales 26]. 



A no-go theorem was formulated by S. Weinberg in ref. [13| which states that scalar 
field cannot naturally solve the vacuum energy problem because one has to impose two 
independent conditions on the adjustment of the potential: U{4>o) — Pvac and the vanishing 
of the derivative of the potential U'(4>o) = at the same point <fi = 0o- However, this 
theorem can be circumvented by more exotic coupling to curvature or by higher spin fields. 



A model with compensating vector field, V^, was proposed in ref. 24J where the La- 
grangian was taken in the form: 

Li = V \F^F^/A + {V.^f} + iRm 2 In f 1 + ^ | . (42) 
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In this theory the time component of the field, V t is unstable in the De Sitter background 
and the solution behaves as: 

V t ~ t + c/t . (43) 

The dominant part of the energy-momentum tensor of this solution is proportional to the 
metric tensor which compensates the original vacuum energy: 

Tp V {V t ) = -pvacg^u + ST^ , (44) 

where 5T^ U tends to zero at large time. 

In this model the gravitational coupling changes with time only logarithmically which 
may possibly agree with the observed bounds on allowed time variation of Gn, but the 
cosmological expansion rate is not connected with the matter content of the universe as is 
dictated by eq. (EJ). 



An interesting possibility opens a massless second rank tensor field S^ u [25] with only 
kinetic term in the action: 

C 2 = ViS a ^S a ^ + V2 S« a S^ + VsS^Sf . (45) 

So it is a free field having minimal gravitational interaction. 

The components of this field in FRW metric satisfy the following equation of motion 

(d 2 + 3Hd t - 6H 2 )S tt - 2H 2 Sjj = , (46) 
(d 2 + 3Hd t -QH 2 )s tJ = 0, (47) 

(d 2 + 3Hdt - 2H 2 )s l3 - 2H 2 5 lJ S tt = , (48) 

where s t j = S t j/a(t) and = Sij/a 2 (t). It is easy to see that time-time component, S tt , 
and isotropic part of space-space component, SV, ~ Sy, are unstable and rise with time in 
De Sitter background. 

Though it is not immediately evident that gravitational coupling to matter in this model 
changes with time, it was argued in ref. that it rises as power of time. 

At this stage the following comment is in order. Normally in a higher spin field theory 
the conditions are imposed eliminating lower spin components, as e.g. in vector field theory 
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the 3- dimensional scalar part (longitudinal component) is eliminated; in second rank tensor 
field the scalar and vector parts are eliminated. Here we have the opposite situation when 
the scalar component of tensor field is physical and probably an elimination of higher spin 
components is necessary. It is unclear if such a theory can be consistently formulated. An 
example of such theory when gauge scalar field is described by t-component of vector is 
presented in paper "Photon and Notoph" by Ogievetsky and Polubarinov j^ . 



A model based on scalar field with "crazy" coupling to gravity was suggested in paper 
and studied in some detail in ref. 



29 1 



301 ] . It is based on the action: 



.4 



\{R + 2A) + Fl {R) + - Ufa R) 



where the system of units with m 2 Pl /%ir = 1 is used. 



The equation of motion for $ has the form: 



D, 



R 



+ u'fa = o. 



(49) 



It is essential that the equation has large coefficient in front of the highest derivative term, 
while quite many modified theories have small coefficient in front of the highest derivative 
and the latter leads to serious problems. 

It is easy to check that solutions of equation (j4"9l asymptotically tends to the state in 
which 



R ~ Pvac + U{<f>) = . 



(50) 



In other words the vacuum energy is eliminated, avoiding the Weinberg no-go theorem. The 
solution has some other nice ("almost realistic") features, for example H = l/2t but this 
result does not depend upon the cosmological matter content and is unstable with respect 
to small fluctuations. To see that one needs to explore the gravitational equation of motion, 
in particular the one for the scalar curvature R, which for a particular choice Fi = C\R 2 has 
the form: 



R + 3{^)\d c 



0) 2 - 4 [Ufa + Pvac ] - 6D 2 



rpn 



(51) 



A desperate attempt to improve the model with non-analytic dependence on R: 



R 2 



R\R\ 



(52) 
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was unsuccessful, though some interesting solutions became stable. 
More general action with a scalar field: 



A — I d 4 xy/—g 



m pi 



(R + 2A) + F 2 (0, R)D lx (j)D' J, (j) + F 3 (0, R)D^D»R - [7(0, R) 



(53) 



167T 

has not yet been explored. Moreover a dependence on R^ v and R^ap can be also included, 
though such terms could lead to ghosts and tachyons in gravitational theory. Anyhow, it is 
difficult to construct a sensible theory without a guiding principle. 



Recently in a series of papers a model with two vector fields has been studied [3l|], which 
according to the author's statement is more successful in the solution of the vacuum energy 
problem then the previous ones. 

Still, we do not have a realistic and consistent with all observations theory explaining 
compensation of the vacuum energy down to the observed small value. As a poor man 
substitution for that, phenomenological approaches are developed. There are two main 
roads to follow: either to introduce into the theory a new field whose equation of state has 
w < —1/3 or to modify gravity at cosmo logically large scales in such a way that it leads to 
accelerated expansion notwithstanding matter. 

V. MODIFIED GRAVITY 

The usual Einstein-Hilbert action is linear in the curvature scalar R. This is the reason 
why the GR equations (T5]) contain, as it is usual in other field theories, only second derivatives 
of metric despite the fact that the action also contains second derivatives. Normally equation 
of motion have one order higher derivatives than the corresponding action, since the kinetic 
term in the latter is a function of the bilinear combination the first order field derivatives. 
Radiative corrections to the GR action lead to generation of higher powers of R in the 
effective action, R n , as well as some invariant combinations of the curvature tensors, e.g. 
powers of R^R^ . This result is true in the limit of small curvature, R <^ m, where m is the 
particle mass in the matter loops in curved space-time. In the opposite limit the effective 



action can be expanded in inverse powers of R. For details one may address book 32]. 

If the action differs from a simple linear GR form, the equation of motion would be 
higher than the second order one. Such equations should contain some pathological features 
as existence of tachyonic solutions or ghosts. However, the theories whose action depends 
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only on a function of the curvature scalar, F(R), are free of such pathologies because, as 
is known, they are equivalent to an addition of a scalar degree of freedom to the usual GR 
with the scalar field satisfying normal second order field equation. That's why modification 
of gravity at large distances is mostly confined to F(R) theories: 

s = ltl d ' x ^ R + F{ M + 5 - • ( 54 ) 



In ref. 



331 ] a model with power law corrections to the Einstein action of the kind 



F(R) = C\R? + R 3 /m| was considered. With m 3 of the order of the neutrino mass, ~ m„, 
such model describes dark energy with the energy density of the order of m 4 , which is quite 
close to the observed value. 

Soon after discovery of the accelerated expansion F(R) theories became quite popular. 
The pioneering works in this direction was done in ref. 34|, which was closely followed by 
ref. 35]. In these works the singular in R action: 

F{R) = -n A /R, (55) 

has been explored with constant parameter /i chosen as fi 2 ~ R c ~ l/tfj to describe the 
observed cosmological acceleration. 

The corresponding equation of motion reads: 

G? + i?) " f G? " if) ^ ~ D ^ GO ^ D - D " (s0 = ■ (56) 

Taking trace over fj, and v of this equation we obtain: 

D 2 R 2 (D a R)(D a R) _R 2 R 4 T R 3 

R 2 ' 1 ' 

Here T = 8irT^ /m 2 pl > 0. This equation has an evident solution in absence of matter R 2 = 

3/i 4 which describes the accelerated De Sitter universe with a constant curvature scalar. 

So far, so good but the small coefficient, /i 4 , in front of the highest derivative or, what 

is the same, the large coefficient, l//i 4 , in front of the non-derivative terms in the presented 

above form of the equation leads to a strong instability in presence of matter 36] . Indeed, 

let us apply equation (j57p for perturbative calculations of the gravitational field of a celestial 

body. 

Let us look for the solution perturbatively with R = Rq + Ri, where Rq is the usual 
solution in the non- modified General Relativity (GR), Rq = —T. Since the gravitational 
field is weak, the flat background metric is assumed. 
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In vacuum, outside the matter source R exponentially tends to zero, i.e to the GR value, 
if /x 4 > 0. So these modified gravity theories agree with the Newtonian limit of the standard 
gravity for sufficiently large \x. 

Now let us consider the internal solution with time dependent matter density. The first 
order correction to the GR curvature, R\ satisfies the equation: 



R 1 -AR 1 -^R 1 + ^f- d J R 1 +R 1 



where {d a Tf = f 2 — (fyTf 



T + 3 



{d a Tf T 3 
T 2 6/x 4 



T 2 

AT H 

2 



T 



(58) 



The last term in the square brackets in the r.h.s. leads to exponential instability of small 
fluctuations as well as to instability of gravitational field created by a regularly varying with 
time mass density of the considered body. The characteristic time of instability is: 



T 



T 3 / 2 



10 



-26 



sec 



-3/2 



/ a 

v g/cm 



(59) 



where p m is the mass density of the body and pT 1 ~ t u ~ 3 • 10 17 sec. This is the dominant 
term in the equation, since T ~ (10 3 sec) -2 (p m /g cm -3 ) and hence the ratio T 3 /p 4 is huge. 

Usually spatial inhomogeneities prevent or inhibit instability, as e.g. happens in the 
case of the Jeans instability. However in our case the term ARi is by far smaller than the 
"unstable" term. The instability would be damped only at the scales of the order of or 
below the Compton wave length of proton. 

To avoid the problem of such instability a modification of the modified grav ity has been 
suggested. We will consider here some class of the models discussed in refs. [37| . Some other 



forms of g ravity modification are reviewed in ref. 38[. The different actions suggested in 



works 



37| have the form: 



Fns(R) 

Fab(R) 
F(R) S 



R v 



\ Rv&c ) 



2n 



- log 

2 6 



XRn 



2 l + 

\ Avac / 

cosh (f - 6) 



2n i 



cosh b 



R 
~2 



(60) 

(61) 
(62) 



Despite different forms these actions result in quite similar consequences. Below we essen- 
tially follow the analysis made in ref. 
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Introducing notation / = R + F(R) we can write the field equations in the form: 

f% - ^; + - D ^)f = || > 

and correspondingly their trace is: 



(63) 



3D 2 f (R) + Rf'(R) - 2f(R) = M^T* . 

The condition of accelerated expansion in absence of matter is that the equation 

Rf'(R) - 2f(R) = 

has the solution R — R± > 0, where R\ is (approximately) constant. 

The following necessary conditions to avoid pathologies are to be satisfied: 

1. Future stability of cosmo logical solutions: 

F\R l )/F"{R l ) > R x . 

2. Classical and quantum stability (gravitational attraction and absence of ghosts): 



(64) 



(65) 



(66) 



F'(R) > 0, F"(R) > 0. 



(67) 



3. Absence of matter instability: 



F'(R) > 0, F"(R) > 0. 



(68) 



Despite considerable improvement these double modified versions still have some serious 
problems. First of all they possess the so called past singularity: in cosmological background 
with decreasing energy density the system must evolve from a singular state with an infinite 
R. In other words, if we travel backward in time from a normal cosmological state, we come 
to a singular state with infinite curvature, while the energy density remains finite. 



Moreover in t 
ity in the future 



re systems with rising mass/energy density the system evolves to a singular- 



4=11 ] . Infinite value of R would be reached in finite (short) time. Following 



ref. 4jJ let us consider version ([62]) in the limit R^> Rq, when one can approximately take: 



F(R) w -XR 



.-if 



(69) 



We analyze the evolution of R in a massive object with time varying density, p ^> p c 
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Gravitational field of such objects is supposed to be weak, so the background metric is 
approximately flat and covariant derivatives can be replaced by the flat ones. Hence equation 
(EH) takes the form: 



{d 2 - A)R - (2n + 2) 



R 2 - (VRf 
R 



+ 



R 2 



3n(2n + 1 



R 



2n 



p2n 
It 



(n + 1) 



2n+l (R + T) = 



(70) 



6n(2n + 1)\R& 

The equation is very much simplified if we choose another unknown function: 
w = F' — —2n\ (Rq/ ' R) 2n+1 which satisfies: 



(d 2 - A)w + U'(w) = 



(71) 



Here potential U(w) is equal to: 

U(w) = \{T- 2XR ) w + ^ 



'7 w 2nu + (u + 



2X\ W 



l+2nu 



2nu 



Q 



2nu 



1 + 2nv 



(72) 



where v = l/(2n + 1) and q = 2n\. 

Notice that the infinite R singularity corresponds to w = 0. 

If only the dominant terms are retained and if the space derivatives are neglected, equation 
[I]) simplifies to: 



™ + T/3-^-^ = 0. 
1 3w v 



(73) 



Potential U would depend upon time, if the mass density of the object changes with time. 
We parametrize it as: 



T = T(t) = T (l + kt) , 

where r is dimensionless time introduced below. 

With dimensionless quantities t — 7T and w = (3z, where 

3q ( i? \ 2(n+1) 



(74) 
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p N 2n+l 



(75) 



the equation further simplifies: 



z" - z~ v + (1 + kt) = . 



(76) 
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One can easily see that the minimum of this time dependent potential moves with time closer 
to z = and becomes more and more shallow. So it is practically evident that independently 
on the initial condition in a system with rising mass density, z reaches zero in finite time and 
correspondingly R = oo. Numerical calculations of ref. 4lj show exactly such a picture. To 
cure this ill-behavior one may add into the action quadratic in curvature term R 2 /6m 2 (42]. 
which does not allow R to reach infinity but the system is stabilized at quite large values of 
R. In the homogeneous case and in the limit of big ratio R/Rq the equation of motion with 
R 2 addition is modified as 



1 - 



R 2n+2 



QXn(2n + 1)R$ 



2n+l m2 



«_(2„ + 2)^- « 2 " + W> 







(77) 



R QXn(2n + l)R 2 n+1 
An important effect which in not taken into account in this equation and which also inhibits 
unbounded rise of R is the particle production by oscillating curvature R. The technique 
for calculations of particles production applicable to the case of modified gravity ( 16011621) 
was worked out in ref. (43] for the case of R 2 gravity in cosmo logical situation (in the 
early universe) where the classical results 44j for particle production were reproduced. The 
developed method is applied to particle production in the modern astronomical systems 
which under certain conditions may be observable sources of cosmic rays 45 1. 

To conclude, there are quite many phenomenological models but no understanding of the 
cosmological acceleration and of the vacuum energy compensation mechanisms are found at 
a deeper level. A very important for future development would be an accurate measurement 
of w. For vacuum energy w must be strictly equal to (-1). The modern data indicate that 
w is close to this value but some deviation, which could be crucial, is possible. 
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